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Abstract. In the framework of ab initio dynamical mean field theory for realistic electronic structure calcu- 
lations a new perturbation scheme which combine the T-matrix and fluctuating exchange approximations 
has been proposed. This method is less computationally expensive than numerically exact quantum Monte 
Carlo technics and give an adequate description of the electronic structure and exchange interactions 
for magnetic metals. We discuss a simple expression for the exchange interactions corresponding to the 
neglecting of the vertex corrections which becomes exact for the spin-wave stiffness in the local approx- 
imation. Electronic structure, correlation effects and exchange interactions for ferromagnetic nickel have 
been discussed. 

PACS. 71.10.-w Theories and models of many-electron systems - 71.15.-m Methods of electronic structure 
calculations 

1 Introduction tion (GGA), which is the base of modern microscopic the- 

ory of solids, is faced with a series of difficulties when 
Electronic structure and magnetic properties of iron-group describing the photoemission, thermoemission and other 
metals are a subject of great interest for a very long pe- spectra of Fe and Ni as well as their finite-temperature 
riod (for review ofearly theories see mi). Density func- magnetic properties (see |,|||,@,| and Refs. therein), 
tional (DF) theory in the form of local spin density ap- The electron correlation effects should be taken into ac- 
proximation (LSDA) or generalized gradient approxima- count to solve these problems. There were a lot of attempts 
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to include these effects in band structure calculations of hole renormalization into the T-matrix which is opposite 
transition metals |l^,|l],|l||l| ,0. Probably the most in this sense to the approach of Ref. 0). 



accurate and successful approach is the use of the dynam- The latter ( "two-step-FLEX" ) approximation has high 
ical mean-filed theory (DMFT, with the solution enough accuracy for the Hubbard model as well as 



of effective impurity problem by numerically exact QMC for real systems with moderate correlations U < W/2 
methods; it has been applied to the magnetism of transi- where U is the Hubbard on-site repulsion energy and W is 
tion metals in Refs.|l^,^. Unfortunately, this approach is the bandwidth [Q. The replacement of the bare Coulomb 
very cumbersome and expensive computationally; besides interaction by the T-matrix can be justified accurately, 
that, the QMC method deals with the "truncated" two- at least for the spin-wave temperature region, both for 
indices interaction matrix instead of the complete four- the Hubbard model and for the s-d exchange (spin- 
indices one (see Q). Alternatively, a scheme has been fermion) model ||2^. However, specific form of the approx- 
proposed in Ref. basing on a multiband spin-polarized imation used in can be improved further by taking into 
generalization of the "fluctuating exchange" (FLEX) ap- account the spin-dependence of the T-matrix. Here we 
proximation by Bickers and Scalapino Original for- present the formulation of the spin-polarized T-matrix- 
mulation of the FLEX approximation treats in the equal FLEX (SPTF) approximation and its application to the 
way both particle-hole (PH) and particle-particle (PP) electron structure of ferromagnetic nickel. 



channels. However, their roles in magnetism are completely 
different: the interaction of electrons with spin fluctuations 
in PH channel leads to the most relevant correlation effects 

0] whereas PP processes are important for the renormal- We start with the general many-body Hamiltonian in the 



2 Computational approach 



AAV 



izations of the effective interactions in spirit of T-matrix LDA+U scheme ||2j 
approach ( "ladder approximation" ) by Galitskii p9| and 
Kanamori . Therefore we used in Ref. a "two-step" 
procedure when, at first, bare matrix vertex is replaced 
by T-matrix, and, secondly, PH channel processes with = - ^ {^i^2\v\ ^1^2) '^ticr'^\2cr'^K'^"^K'^ ' 

{Ai}crcr' 

where A — im are the site number (i) and orbital (m) 
quantum numbers, a , i is the spin projection, c^,c 
are the Fermion creation and annihilation operators, Ht 
is the effective single-particle Hamiltonian from the LDA, 
corrected for the double-counting of average interactions 



this effective interaction are taken into account explicitly. 
Note that the first attempt to combine the T-matrix and 
particle-hole correlations in relation with the problem of 
the magnetism of transition metals has been done by Lieb- 
sch B ) but in a different way (introducing of the particle- 



(12 \v\ 34) = J drdrVi*(r)^2*(r')« (r " r')Mr)Mr^'), 

(2) 

where we define for briefness Ai = 1 etc. Following Ref. 
]l9| we take into account the ladder (T-matrix) renormal- 
ization of the effective interaction: 
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among correlated electrons as it will be described below, functions of a ferromagnet, at least, for spin-wave tem- 
and the Coulomb matrix elements are defined in the stan- perature region ||2^,^). We improve the approximation 
dard way by taking into account the T-matrix spin dependence. 

When considering the particle-hole channel we replace in 
the Hamiltonian (0) u ^ T°'°' which is the solution of 
Eq.(||) ai f2 — 0. Eq. (Q) is exact in the limit of low elec- 
tron (or hole) density which is important for the criterion 
of magnetism, e.g., in the case of nickel (with almost com- 
pletely filled d- band). 

13 T--' izf2) 24) = (13 \v\ 24) - ^ J] J] (13 \v\ 57) * 

f uj 5678 Now we rewrite the effective Hamiltonian (nf) with the 

Gl, {^.) G?; {^n - ^.) (68 \t^^' m\ 24) , (3) ,,pi„,, |,| 34^ ^y (l2 |t-' | 34) in Hu- To con- 

where w = (2n + l)7rT are the Matsubara frequencies for sider the correlation effects due to PH channel we have 

temperature T = /J^^ (n = 0, ±1, ...). Further we rewrite separate density [d) and magnetic (m) channels as in 

the perturbation theory in terms of this effective interac- ^g^- fl^ 
tion matrix. 

At first, we take into account the "Hartree" and "Fock" 

diagrams with the replacement of the bare interaction by '^12 = ^c^|C2| -I- c^C2x 

1 



the T-matrix = _ (^c+ cat - c+ C2i 

^{2,^ (*^) = 4 E E (13 \t""' m\ 24) - ""^i ^ '^T^2i 

^ a 34a' 

^iZ' (-) = E E \T^^ (mm Gl, {^n - ^.) ^ '^^''^ ' 

^ n 34 

(4) 

Then the interaction Hamiltonian can be rewritten in the 
Note that 17^"^^-' -I- E'^'^^^ contains exactly all the second- following matrix form 
order contributions as it can be easily seen from the cor- 
responding Feynman diagrams. Now we have to consider ~ 2"^' * V^" * D -|- m * * m -|- w * * n 
the contribution of particle-hole excitations to the self- 

(6) 

energy Similar to we will replace in the corresponding ^^ere * means the matrix multiplication with respect to 
diagrams the bare interaction by the static hmit of the T- P^i'^^ °^ orbital indices, e.g. 

matrix (as it was already mentioned, it can be justified by 

the explicit calculation of the electron and magnon Green ■ ^ 
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the supervector D defined as 



D = {d,m°) ,D+ = 



and the effective interactions have the foUowing form: 

(K;^)n^2. = -(l2|T^^|2'l') 

^dd ^dm 
^md ^dd 

^,22' = ^ E \t"^' I 1'2') - ^ E 2'!') 

n'> = J E T^"' 1'2') - I E 2'!') 



(11) 



Y 



and the matrix of bare longitudinal susceptibility is 

in the dd-, dmP-, mPd-, and mPmP- channels (d, m° = 
1,2 in the supermatrix indices). An important feature of 
these equations is the coupling of longitudinal magnetic 
fluctuations and of density fluctuations. It is not present 
in the one-band Hubbard model due to the absence of the 
interaction of electrons with parallel spins. For this case 



^11'. 22' — ^22M1' 



md 



1 



Eqs. 10) coincides with the well-known result ||25| . 

Now we can write the particle-hole contribution to the 
self-energy. Similar to Ref.H one has 



-[(12 |rTT| i'2') _ (12 |T^^| 1'2') - (12 \T^^\ 1'2') + 
(12 |T^T| i'2') _ (12 |tTT| 2'1') + (12 |T^^| 2'1')] 



(7) 



(12) 



34. cr' 



with the P-H fluctuation potential matrix: 



To calculate the particle-hole (P-H) contribution to the 
electron self-energy we first have to write the expressions 
for the generalized susceptibilities, both transverse a-nd 
longitudinal x". The corresponding expressions are the 
same as in Rcf . [Q but with another definition of the inter- 
action vertices. One has 

x+- {lu) = [1 + v„\ * r^^iito)] {ico) , (8) 

where 



(13) 



VFTT (iuj) (icu) 
(iuj) W^^ (iuj) 

where the spin-dependent effective potentials are defined 



as 



\Yll ^ _y\\ ^ 
2 



x" - xl! 



(14) 



(t) = -G^3(^)G4^i i-r) 



(9) 



is an "empty loop" susceptibility and r{iLu) is its Fourier 
transform, r is the imaginary time. The corresponding 
longitudinal susceptibility matrix has a more complicated 
form: 

x" (tiu) =\l + * x|l(*^)l ' * xll(«^), (10) 



where x" , Xo differ from x" , Xg ^^'^ replacement of .T^^ <^ 
pll in Eq. (prj) . We have subtracted the second-order con- 
tributions since they have already been taken into account 
inEq.d). 

Our final expression for the self energy is 

E = S^^"^ + S^^^^ + S^P"^ (15) 
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This formulation takes into account accurately spin-polaron The Weiss field function Q is required to obey the self 
effects because of the interaction with magnetic fluctua- consistency condition, which restores translational invari- 



tions 1 22, 2^, the energy dependence of T-matrix which ance to the impurity model description: 

is important for describing the satellite effects in Ni B, 

" G,M = ^[(ztu + M)l-iJ(k)-SfM]-i (18) 

contains exact second-order terms in v and is rigorous (be- k 

cause of the first term) for almost filled or almost empty ^'^^^^ ^ chemical potential, H{k) is the LDA Hamil- 

^g^^jg tonian in an orthogonal basis. The local matrix 'Ef^ is the 

mi T^T T-i-sr • 1 1 • • • 1 T sum of two terms, the impurity self energy and a so-called 

1 he i" LhjX approximation can be used m principle di- t j oj 

, , , , , , , , , , , . . , , , , "double counting " correction, Erir which is meant to sub- 

rectly to the crystal problem taking into account the mo- 

, , „ , , , „ , . , , , , , tract the average electron-electron interactions already in- 

mentum dependence of the self-energy, which would lead 



eluded in the LDA Hamiltonian. For metallic systems we 



to very cumbersome calculations. To overcome this com- 

... 1 • n -c iRi 1 1 propose the general form of dc-correction: S'^'^ (zcj) = So- (iix))- 

putational problem, we use as m Kef.|[7| a local approx- ^ ^ o a \ / o \ , 

^. ^ ,„ 1- i 1 • lTrc'S„ (0) for non-magnetic LDA Hamiltonian S and 

imation to the self-energy, corresponding to combination 2 " " / => o 

of the SPTF approach presented above with the DMFT ^'^^ = ^'^^'^^ LSDA Hamil- 



tonian. This is motivated by the fact that the static part 
of the correlation effects are already well described in 
the density functional theory. Only the d-part of the self- 



theory. The latter maps the many-body system onto a 
multi-orbital quantum impurity, i.e. a set of local degrees 

of freedom in a bath described by the Weiss field function 

^ mi • /1 ^ \ r ^ M ■ J- energy is presented in our calculations, therefore Jl^'^ = 

y. ihe impurity action (here c(t) = [Cma(T)\ is a vector oj > o- 

r ^ • ui \ • • u for s- and p- states as well as for non-diagonal d — s.p 

of Grassman variables) is given by: o jj- 

contributions. 

r/3 /-/J 



/ / 1 

^"^^ ^ Jo ^'^ Jo ^'^'^^^^ ^'^^^ ('r>^')c(T-')] + In spirit of the DMFT approach we have to use the 

/ drHu [c+(t),c(t)] (16) Weiss function instead of Ga in all the expressions 

Jo 

when calculating the self-energy on a given site. Similar to 

It describes the spin, orbital, energy and temperature de- 

the one-band DMFT-perturbation scheme [^,[2^ we keep 

pendent interactions of particular magnetic 3o?-atom with 

the static mean-field term in the bath Green functions: 

the rest of the crystal and is used to compute the local 

g-\iuj) - G-\iuj) + S,(zc^) - (0). 

Greens function matrix: 

G„{t - t') = -^J D[c,c+]e-^^ffc{T)c+{T') (17) 3 Electronic structure of nickel 

{Z is the partition function) and the impurity self energy We have started from the non spin-polarized LDA or spin- 
Q~^{iuj) — G~^{iuj) — ■ polarized LSDA band structure of nickel within the TB- 
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LMTO method |^ in the minimal s,p, d basis set and density of states (DOS) on U. To find DOS we apphed a 

used numerical orthogonalization to find the H(k.) Hamil- Pade approximant method pTf for the analytical continua- 

tonian in Eq. (^8|) . We take into account of the Coulomb tion of the Green function from the Matsubara frequencies 

interactions only between d-states. The correct parame- to the real energy axis. To find the self-consistent solution 

terization of the i/j/part is indeed a serious problem. For of the SPTF equations we used 1024 Matsubara frequen- 

example, first-principle estimations of average Coulomb cies and the FFT-scheme with the energy cut-off at 25 eV 



interactions ([/) [p9|, 13 lead to unreasonably large value and temperature around 200 K. The sum over irreducible 
of order of 5-6 eV in comparison with experimental val- Brillouin zone have been made with 256 k-points. 
ues of the [/-parameter in the range of 1-2 eV for iron||l^. 
Semiempirical analysis of the appropriate interaction value 

Comparison of the LDA density of state and the SPTF 

pot gives [/ ~ 3 eV. It is shown in Refs. that an ad- 
calculation with DMFT self-consistency for the local self- 
equate description of a broad circle of the properties of 

energy matrix (Fig.l) shows that the latter does repro- 

Fe and Ni in the LDA+DMFT scheme is possible when 

duce the three most important characteristic features of 

choosing [/ ~ 2 — 3 eV. The difficulties with an ab initio 

correlation effects for nickel: 6 eV satellite, 30% narrowing 

determination of the correct value of U are connected with 

of the d-bandwidth and 50% reduction of exchange split- 
complicated screening problems, definitions of orthogonal 

tings in comparison with the LSDA band structure 

orbitals in the crystal, and contributions of the intersite in- 

^,|^|3^. For U—2 eV the position of satellite is repro- 

teractions. In the quasiatomic (spherical) approximation 

duced quite well, while for U—3 eV it is shifted to the 

the full [/-matrix for the d— shell is determined by the 

lower energies. Note that the LDA-I-DMFT consideration 

three parameters U, J and SJ or equivalently by effective 

with the QMC solution of the effective impurity problem 

Slater integrals F^, F^ and |||. For example, U= F°, 

gives an adequate description of the electronic structure 

J=(F^+F'*)/14 and we use the simplest way of estimating 

of Ni for the choice U — 3 eV ||] . The narrowing of the 

SJ or F* keeping the ratio F^/F* equal to its atomic value 

d-bandwidth in our calculations is reasonable for the both 

0.625 0. 

[/-values. The non-magnetic LDA starting Hamiltonian is 
Note that the value of intra-atomic (Hund) exchange better than the LSDA one for correct description of the 
interaction J is not sensitive to the screening and approxi- 50% reduction of the spin-splittings in nickel, while for 
mately equals 1 eV in different estimations [p9| ; further we magnetic LSDA Hamiltonian the the spin-splitting in the 
have chosen J — 1 eV. For the most important parameter quasiparticle DOS remains approximately the same like 
[/, which defines the bare vertex matrix (Eq. (||)), we took in the LSDA results (Fig.l). The local magnetic moment 
the values U = 2 and 3 eV to check the dependence of the on nickel atom is not very sensitive to the U-values and 
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account in our scheme due to presence of the off-diagonal 
terms in the effective potential Eq.(0|). 




Energy, eV 



Fig. 1. Spin- up (full lines) and spin-down (dashed lines) den- 
sity of d-states for ferromagnetic nickel in the LSDA and the 
LDA+SPTF (LSDA+SPTF) calculations for different aver- 
age Coulomb interaction U with J = 1 eV and temperature 
T=200 K. 

is equal to 0.56 fiB for U=2 eV LDA+SPTF and 0.58 /is 
for U=3 eV LSDA-hSPTF calculations. 

Another important correlation effect is an essential re- 
duction of the spin polarization near the Fermi level in 
comparison with the LSDA calculations. This is connected 
with the spin-polaron effects because of the mixing of the 
spin- up and spin-down states p6| . They are taken into 



4 Exchange interactions in nickel 

Calculating the variation of the thermodynamic poten- 
tial with respect to small spin rotations with the use of 
the "local force theorem" an effective exchange interaction 
parameters can be found in the following form [T^ 

J,, = -rr^L(^f4i:|4) (19) 

where Z'f ~ ^ (^sj — sj^^ . Correspondingly, the magnon 
dispersion relation Wq for a ferromagnet is defined by the 
formula 

= ^ [J(0) - J(q)] (20) 

where M is the magnetic moment per unit cell, J(q) is 
the Fourier transform of the exchange integrals defined by 
Eq.(^9|). The expression for the stiffness tensor Da^, 

= Dafiqaqp, q ^ 0, (21) 

reads 

^.. = -^T..,E(^^^^^^) (22) 

These results generalize the LSDA expressions of Rcf. 
to the case of correlated systems. One can show (see Ap- 
pendix) that they can be derived using a standard dia- 
gram approach under two assumptions: (i) the locality of 
the self-energy (which is fulfilled in the DMFT) and 
(ii) the neglecting of the vertex corrections. The expres- 
sion (|2|) for the stiffness constant turns out to be exact 
in the framework of the DMFT. 
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5 Conclusions 

Here we have presented the results of new SPTF approxi- 
mation in the framework of first-principle dynamical mean 
field theory for magnetic transition metals. This approxi- 
mation combining the T-matrix and FLEX schemes gives 
a satisfactory description of both electronic and magnon 
spectra of Ni. In contrast with the QMC method for the 
solution of the effective impurity problem, this approach, 
being less rigorous, is not so time- and resource-consuming 
and allows to work with the most general rotationally in- 
variant form of the Coulomb on-site interaction. 
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Fig. 2. Spin-wave spectrum for ferromagnetic nickel in 
LDA+SPTF scheme with (7 = 2 eV and J = 1 eV in compar- 
ison with experimental magnon spectrum (Ref. 36) in _r — L 
direction. 

We have calculated the magnon spectrum for the op- 
timal choice J7=2 eV and J—\ eV using SPTF calcula- 
tions with the non-magnetic LDA as a starting point. The 
computational results are shown in Fig. 2; the calculated 
spin- wave stiffness constant for Ni is found to be 13 = 450 
meV/A^ in an excellent agreement with the experimental 
value of 455 meV/A'^ Note that simple approximation 
for exchange interactions is not allowed us to investigate 
the problem of the optical mode in magnon spectrum of 
nickel Eol. 
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Appendix: Exchange interactions and vertex 
corrections 

In order to elucidate the approximation behind the ex- 
pression for the exchange parameters (Eq. p^ , we con- 
sider the energy of a spiral magnetic configuration with 
the rigid rotation of the spinor-electron operators on site 
i by the polar angles Qi and (/s^: 

Cim ^ U (^ii ^i) ^im 

where 

cos 0/2 sin6'/2exp (— i(/3) 
sin 6*72 exp (i(/3) cos6'/2 
assuming that 6i = const and tpi — qRi where Ri is the 
site lattice vector. Since we take into account only on-site 



u{e,ip) = 
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correlation effects tlie interaction term in the Hamiltonian vertex are local. This can be proven, e.g., using the Ward- 
is invariant under that transformation, and the change of Takahashi identities ||3^ . Therefore, the expression for the 
the Hamiltonian is stiffness constant of the ferromagnet (Eq. (p2|)) appears to 
SH TvLa [Uj4 (U^Uj - 1) Cj] = SiH + 52H be exact in the framework of DMFT |3|. 

SiH = sin2 f ^ TvLa [{t (k+q) - t (k)) c+c^] 

k 

52H = ^sinej^TrL Ujcf^CjT^ (exp (iqRj) - exp (iqRj)) References 

ij 

(23) 

Consider further the case of small 0, we can calculate ^- ^- Herring Magnetism, vol. 4 (Academic Press, New York, 
the variation of the total energy to lowest order in 9 which '^'^'^ 

corresponds to the first order in S,H and the second order 2- S. V. Vonsovsky, Magnetism (Wiley, New York, 1974). 



in 69H 



3. T. Moriya, Spin Fluctuations in Itinerant Electron Mag- 

2 netism (Springer, Berlin, 1985). 

6E - -Tr^ V[i(k + q) - t(k)]{n^ + 

4 ^ 4. V.Yu.Irkhin, M.I. Katsnelson, and A. V.Trefilov, J. Phys.: 

Tr^lik, q)Gdk + ,)[t(k + q) - t(k)]GT (fc)]}, (24) ^^^^^^^ ^^^^^^ ^ ^^^3 ^^^^^^^ 

where rik = TrL^(c^Ck), q, k are four- vectors with com- 5. D. VoUhardt, N. Bluemer, K. Held, M. KoUar, J. Schlipf, M. 

ponent (q, 0) and (k,icj), 7 is the three-leg vertex. Our Ulmke, and J. Wahle, in: Advances In Solid State Physics, 

main approximation is to neglect the vertex corrections P' (Vieweg, Wiesbaden 1999). 

(7 = 1). In this case the previous equation takes the fol- Band-Ferromagnetism. Ground State and Fimte- 

, . J. Temperature Phenomena , ed. by K. Baberschke, M. 

lowing form: 

2 Donath, and W. Nolting (Springer, Berlin, 2001). 

SE=-TrL^Y.^t{k + q)-t{k)]* (25) 

4 7. M. I. Katsnelson and A. I. Lichtenstein, J. Phys.: Cond. 

G,{k + q)[Gl\k + q)- G^\k) + t{k + q) - m]G^(k) ^^^^^^ ^^3^ ^^^^^^ 

Using the following consequence of the Dyson equation: 8. A. I. Lichtenstein, M. I. Katsnelson, and G. Kotliar, Phys. 

Rev. Lett. 87, 067205 (2001). 

t (k+q) - t (k) = Gf 1 (fc) - {k + q) + {E) - [E] 

9. A. Liebsch, Phys. Rev. B 23 , 5203 (1981). 
one can rewrite the Eq.(|2|) in the form: 5E = ^[J(0) - ^ ^^^g^^^ p Ducastelle, and D. Spanjaard, J. Phys. 

J(q)] with the exchange integrals corresponding to Eq. (Paris) 43 341 (1982) 

(|l|). We conclude that the expression for J„ is accu- p Manghi, V. Bellini, and C. Arcangelli, Phys. Rev. B 

rate if the vertex corrections can be neglected. Note that 55, 7149 (1997); F. Manghi, V. Bellini, J. Osterwalder, T. J. 
the hmit of smaU q this neglecting can be justified rigor- Kreutz, P. Aebi, and C. Arcangeh Phys. Rev. B 59, R10409 
ously, provided that the self-energy and three-leg scalar (1999). 



10 M. I. Katsnelson, A. I. Lichtenstein: Electronic structure and magnetic properties of correlated metals 

12. W. Nolting, S. Rex, and S. Mathi Jaya, J. Phys. : Condens. 27. Y. H. Szczech, M. A. Tusch, and D. E. Logan, Phys. Rev. 
Matter 9, 1301 (1987). Lett. 74, 2804 (1995). 

13. M. M. Steiner, R. C. Albers, and L. J. Sham, Phys. Rev. 28. O. K. Andersen and O. Jepsen, Phys. Rev. Lett. 53, 2571 
B 45, 13272 (1992). (1984). 

14. V. Drchal, V. Janis, J. Kudrunovsky, Phys. Rev. B 60, 29. V. I. Anisimov and O. Gunnarsson, Phys. Rev. B 43, 7570 
15664 (1999). (1990). 

15. W. Metzner and D. VoUhardt, Phys. Rev. Lett. 62, 324 30' A- M. Oles and G. StoUhoff, Phys. Rev. B 29, 314 (1984). 
(1989); E. MxiUer-Hartmann, Z. Phys. B 74, 507 (1989); A. H- J- Vidberg and J. W. Serene, J. Low Temp. Phys. 29, 
Georges and G. Kotliar, Phys. Rev. B 45, 6479 (1992); T. ^"^^ (1977). 

Pruschke, D. Cox, and M. Jarrell, Phys. Rev. B 47, 3553 ^2. M. Iwan, F. J. Himpsel, and D. E. Eastman, Phys. Rev. 



(1993). 



Lett. 43, 1829 (1979) 



T„ . „ ^1- ur J AT r. u 33. W. Eberhardt and E. W. Plummer, Phys. Rev. B 21, 3245 

16. A. Georges, G. Kothar, W. Krauth, and M. Rozenberg, 

(1980). 



Rev. Mod. Phys.68, 13 (1996 ). 

17. M. L Katsnelson and A. I. Lichtenstein, Phys. Rev. B 61, 
8906 (2000). 

18. N. E. Bickers and D. J. Scalapino, Ann. Phys. (N.Y.) 193, 
206 (1989). 

19. V. M. Galitskii, Zh. Eksper. Teor. Fiz. 34, 151, 1011 (1958) 

20. J. Kanamori, Prog. Theor. Phys. 30, 275 (1963). 

21. M. Fleck, A. I. Liechtenstein, A. M. Oles, L. Hedin, and 
V. I. Anisimov, Phys. Rev. Lett. 80, 2393 (1998). 

22. V. Yu. Irkhin and M. L Katsnelson, J. Phys. : Condens. 
Matter 2, 7151 (1990). 

23. V. Yu. Irkhin and M. I. Katsnelson, Sov. Phys. - JETP 61, 
306 (1985). 

24. V. I. Anisimov, F. Aryasetiawan, and A. I. Lichtenstein, 
J. Phys.: Condens. Matter 9, 767 (1997). 

25. T. Izuyama, D. Kim, and R. Kubo, J. Phys. Soc. Japan 
18, 1025 (1963). 

26. V. Yu. Irkhin and M. I. Katsnelson, Physics - Uspekhi 
37,659 (1994). 



34. K.-P. Kaemper, W. Schmitt, and G. Guentherodt, Phys. 
Rev. B 42, 10696 (1990). 

35. K. N. Ahmann, D. Y. Petrovykh, G. J. Mankey, N. Shan- 
non, N. Gilman, M. Hochstrasser, R. F. Willis, and F. J. 
Himpsel, Phys. Rev. B 61, 15661 (2000). 

36. A. I. Liechtenstein, M. I. Katsnelson, and V. A. Gubanov, 
J. Phys. F 14, L125 (1984); A. I. Liechtenstein, M. I. Kat- 
snelson, V. P. Antropov, and V. A. Gubanov, J. Magn. Magn. 
Mater. 67, 65 (1987). 

37. H. A. Mook and D. McK. Paul, Phys. Rev. Lett. 54, 227 
(1985) 

38. J. A. Hertz and D. M. Edwards, J. Phys. F 3, 2174 (1973). 

39. A. I. Lichtesntein and M. I. Katsnelson, In: Ref. Q, p. 75. 

40. S.Y. Savrasov, Phys. Rev. Lett. 81, 2570 (1998). 



